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Power  losses  and  voltage  instability  are  major  problems  in  present  power  systems.  It  has  become  more  complex 
day  by  day  due  to  less  security  and  reliability.  Flexible  AC  transmission  systems  (FACTS)  controllers  have  been  mainly 
used  for  solving  various  power  system  steady  state  control  problems.  Flexible  AC  transmission  systems  or  FACTS  are 
devices  which  allow  the  flexible  and  dynamic  control  of  power  systems  and  enhancement  of  system  stability  using  FACTS 
controllers.  Based  on  PSO  algorithm  is  an  effective  method  for  finding  the  optimal  choice  and  location  of  FACTS 
controllers.  It  also  increases  load  ability  of  the  line  and  minimizes  losses.  This  paper  presents  comparative  study  of  GA  and 
PSO  algorithms  for  one  of  the  FACTS  controller  i.e.,  UPFC  device.  The  suggested  algorithm  has  been  applied  to  IEEE-30 
bus  system. 
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The  power  system  stability  is  considered  as  one  of  the  most  significant  concepts  of  power  systems  quality. 
Nonlinearity  characteristics,  high  complexity  and  time  varying  behaviour  of  power  systems  have  considered  widespread 
challenges  to  stability  of  the  power  systems.  There  are  different  problems  in  maintaining  the  stability  of  power  systems. 

Voltage  stability  refer  to  the  ability  of  power  system  to  maintain  steady  voltages  at  all  buses  in  the  system  after 
being  subjected  to  a  disturbance  from  a  given  initial  operating  condition.  Voltage  instability  is  loss  of  load  in  an  area,  or 
tripping  of  transmission  lines  and  other  elements  by  their  protective  systems  leading  to  cascading  outages.  Loss  of 
synchronism  of  some  generators  may  result  from  these  outages  or  from  operating  conditions  that  violate  field  current  limit. 
So  that  control  of  load  flow  in  order  to  have  more  efficient,  reliable,  and  secure  system  is  in  the  interest  of  the  transmission 
system  operator  (TSO).  To  overcome  this  problem  FACTS  device  are  introduced. 

FACTS  can  increase  controllability  and  optimize  the  utilization  of  the  existing  power  system  capacity.  FACTS 
employ  high  speed  thyristors  for  switching  in  or  out  transmission  line  components  such  as  capacitor,  reactor  or  phase 
shifting  transformer  for  some  desirable  performance  of  the  systems.  UPFC  is  the  one  of  the  most  important  device  in 
FACTS  controllers.  It  is  in  multiple  line  compensation  are  integrated  into  a  generalized  power  flow  controller  that  is  able 
to  maintain  prescribed,  and  independently  controllable,  real  and  reactive  power  flow  in  the  line. 

A  Genetic  algorithm  (GA)  is  a  search  heuristic  that  mimics  the  process  of  natural  selection.  This  heuristic 
(also  sometimes  called  metaheuristic)  is  routinely  used  to  generate  useful  solutions  to  optimization  and  search  problems. 

Particle  Swarm  Optimization  (PSO)  is  an  extremely  simple  algorithm  that  seems  to  be  effective  for  optimizing  a 
wide  range  of  functions.  PSO  is  applied  for  solving  various  optimization  problems  in  electrical  engineering. 
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PSO  shared  many  similarities  with  evolutionary  computation  techniques  such  as  Genetic  Algorithms  (GA). 
The  system  is  initialized  with  a  population  of  random  solutions  and  searches  for  optima  by  updating  generations.  However, 
unlike  GA,  PSO  have  no  evolution  operators.  In  PSO,  the  potential  solutions,  called  particles,  fly  through  the  problem 
space  by  following  the  current  optimum  particles.  PSO  is  easy  to  implement  and  there  are  few  parameters  to  adjust. 

MODELLING  OF  UPFC 

Unified  Power  Flow  Controller  is  the  one  of  main  component  of  the  FACTS  devices.  It  consists  of  two  switching 
converters  that  are  operated  from  a  common  direct  current  link  provided  by  direct  current  storage  capacitor.  Converter  2 
provides  the  main  function  of  UPFC  by  injecting  an  alternative  current  voltage  of  controllable  magnitude  and  phase  angle 
in  series  with  the  transmission  line  via  a  series  transformer  1 .  The  basic  function  of  converter  1  is  to  supply  or  absorb  the 
real  power  demand  by  converter  2  at  the  common  direct  current  link.  The  active  current  component  of  converter  1  is 
obtained  from  power  balance  between  the  series  and  the  shunt  converters  and  the  reactive  component  can  be  independently 
controlled  to  provide  necessary  voltage  support  at  the  bus. 
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Figure  1:  Structure  of  UPFC 

Power  flow  through  the  transmission  line  i-j  namely  is  depend  on  line  reactance^- ,  bus  voltage  magnitudes 


Vt,  Vj  and  phase  angle  between  sending  and  receiving  buses  St  —  Sj.  This  is  expressed  by  Eq.  1. 
Py  =  —±sm{8i-8j) 

■A  i  j 
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Two  types  of  UPFC  models  are  available.  One  is  a  coupled  model  and  other  is  decoupled  model.  In  the  first  type, 
UPFC  is  modelled  with  series  combination  of  a  voltage  source  and  impedance  in  the  transmission  line.  In  decoupled 
model,  UPFC  is  modelled  with  between  separated  buses.  First  model  is  more  complex  compared  with  the  second  one 
because  modification  of  Jacobian  matrix  in  coupled  model  is  evitable.  In  decoupled  model  can  be  easily  implemented  in 
conventional  power  low  algorithm  without  modification  of  Jacobian  matrix  elements.  Decoupled  model  has  been  used  for 
modelling  UPFC  in  power  flow  study.  To  obtain  UPFC  model  in  load  low  study,  it  is  represented  by  four  variables 
(i.e.  Pul,  Qul,  Pu2  and  QU2).  Assuming  UPFC  to  be  lossless  and  real  power  flow  from  bus  i  to  bus  j  can  be  expressed  as:  [8] 
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Figure  2:  Decoupled  Model  of  UPFC 
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Although  UPFC  can  control  the  power  flow,  but  cannot  generate  the  real  power.  So 


Pui  +  PU2  =  0 


(3) 


It  is  assumed  that  the  UPFCs  are  installed  in  the  middle  of  lines.  Each  reactive  power  output  of  UPFC  Qul  Qu2  can 
be  set  to  an  arbitrary  value  depends  on  rating  of  UPFC  to  maintain  bus  voltage. 

Optimal  Location  for  UPFC 

In  order  to  find  the  optimal  location  or  the  UPFC  to  be  placed  the  bus  which  is  mostly  affected  during  faults  has  to 
be  identified.  With  the  increased  loading  of  transmission  and  distribution  lines,  voltage  instability  problem  has  become  a 
concern  and  serious  issue  for  power  system  planners  and  operators.  The  main  challenge  of  this  problem  is  to  narrow  down 
the  locations  where  voltage  instability  could  be  initiated  and  to  understand  the  origin  of  the  problem.  One  effective  way  to 
narrow  down  the  workspace  is  to  identify  weak  buses  in  the  systems,  which  are  most  likely  to  face  voltage  collapse  and 
transmission  line  losses. 

VOLTAGE  STABILITY  INDEX  (VSI) 

Consider  the  power  network  where  n  is  the  total  number  of  buses  with  1,  2...g  generator  buses,  and  g+l...n 
remaining  (n-g)  buses.  For  a  given  system  operating  condition,  using  the  load  flow  results,  the  voltage  stability  index  'L' 
can  be  calculated  as  [7] 


Where  j=g+l....n  and  all  the  terms  inside  the  sigma  on  the  right  hand  side  complex  quantities.  The  complex 
values  of  Fj.are  obtained  from  Ybus  matrix  of  power  system.  For  a  given  operating  condition. 


Where  IG,IL  and  VG,VL  represent  complex  current  and  voltage  vectors  at  the  generator  nodes  and  load 
nodes  [yGG],[yGL],  |YLG]and|YLL]  are  corresponding  partitioned  portions  of  the  YBUS  matrix. 

This  analysis  will  be  carried  out  only  for  the  load  buses;  hence  the  index  obtained  is  for  load  buses  only. 
For  stability  the  index  L  must  not  be  more  than  1  for  any  of  the  nodes  j.  the  global  index  for  stability  of  the  given  power 
system  is  defined  to  be  L=  maximum  of  Lj  for  all  j  (load  buses).  The  index  far  away  from  1  and  close  0  indicates  voltage 
stability.  The  L  index  will  give  the  scalar  number  to  each  load  bus.  Among  the  various  indices  for  voltage  stability  and 
voltage  collapse  prediction  (i.e.  far  away  from  1  and  close  to  1  or  >1  respectively),  the  L  index  will  give  more  accurate 
results.  The  L  indices  for  given  load  conditions  are  calculated  for  all  load  buses  and  the  maximum  of  the  L  indices  gives 
the  proximity  of  the  system  to  voltage  collapse. 

EVOLUTION  OF  OPTIMIZATION  TECHNIQUES 
Genetic  Algorithm 

GAs,  as  powerful  and  broadly  applicable  stochastic  search  and  optimization  techniques,  are  perhaps  the  most 
widely  known  types  of  evolutionary  computation  methods  today.  GA  of  evolutionary  computation  is  inspired  by  Darwin's 
theory  of  evolution.  A  GA  generates  a  population  of  possible  solutions  encoded  as  chromosomes,  evaluates  their  fitness 
and  creates  a  new  population  by  applying  genetic  operators  which  are  selection,  crossover  and  mutation. 


Lj  =  i-2: 


(4) 


YGG  YGLUVG] 
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•  Selection 

In  proposed  GA,  method  of  tournament  selection  is  used  for  selection.  This  method  chooses  each  parent  by 
choosing  nt  (tournament  size)  players  randomly  and  choosing  the  best  individual  out  of  that  set  to  be  a  parent.  In  this  paper 
ntis  chosen  nt  =  4. 

•  Cross  Over 

Cross  over  allows  the  genes  from  different  parents  to  be  combined  in  children  by  exchanging  materials  between 
two  parents.  Cross  over  function  randomly  selects  a  gene  at  the  same  coordinate  from  one  of  two  parents  and  assign  it  to 
child.  For  each  chromosome,  a  random  number  is  selected.  If  this  number  is  between  0.01  and  0.3,  two  parents  are 
combined;  else  chromosome  is  transferred  with  no  cross  over. 

•  Mutation 

GA  creates  mutation  children  by  randomly  changing  the  genes  of  individual  parents.  In  this  paper,  GA  adds  a 
random  vector  from  a  Gaussian  distribution  to  the  parents.  For  each  chromosome,  random  number  is  selected.  If  this 
number  is  between  0.01  and  0.1,  mutation  process  is  applied;  else  chromosome  is  transferred  with  no  mutation. 

Particle  Swarm  Optimization 

Inspired  by  social  behaviour  of  bird  flocking  or  fish  schooling,  Eberhart  and  Kennedy  first  developed  the  particle 
swarm  optimization  (PSO)  algorithm  in  1995.  PSO,  as  a  branch  of  evolutionary  computation,  has  been  successfully  applied 
in  many  research  and  application  areas  in  the  past  several  years.  The  PSO  as  an  optimization  tool  provides  a 
population-based  search  procedure  in  which  individuals  called  particles  change  their  position  (state)  with  time.  In  a  PSO 
system,  particles  fly  around  in  a  multidimensional  search  space.  During  flight,  each  particle  adjust  its  position  according  to 
its  own  experience  (This  value  is  called  Pbest),  and  according  to  the  experience  of  a  neighbouring  particle 
(This  value  is  called  Gbest),  made  use  of  the  best  position  encountered  by  itself  and  its  neighbour.  [5]  Each  particle  is 
treated  as  a  point  in  a  D  dimensional  space. 

The  ith  particle  is  represented  as 

%i  —  (Xil>Xi2>  %id) 

The  best  previous  position  of  the  ith  particle  is  recorded  and  represented  as 

Pbesti  ~  (Pbestil,Pbesti2i  ^bestid) 

The  index  of  the  best  particle  among  the  entire  particle  in  the  population  is  represented  by  the  symbol  gbestid- 
The  rate  of  the  position  change  (velocity)  for  the  particle  is  represented  as 

vt  =  (ytl,vt2  vid) 

After  the  search  procedure  for  the  Pbest  id  and  gbest  id  the  modified  velocity  and  position  of  each  particle  can  be 
calculated  using  the  current  velocity  and  the  distance  from  P^estid  to  dbestid  as  shown  in  the  following  equation. 

Vid  (t  +  l)  =  W*  Vid  (t)  +  d  *  randv  (Pbestid  -  Xid  (t))  +  C2  *  rand2  (gbestid  -  Xid  (t))  (5) 
Where 

i=  1,2,  3...n 
d=  1,2,3. ...m 
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n:  Number  of  particles 

m:  Solution's  dimensions  i.e.,  the  number  of  control  variables 

V id(0:  Current  velocity  of  the  i  at  iteration  t 

Vid(t  +  1):  Modified  velocity  of  particle  i 

Rand:  Random  number  between  0  and  1 

Xid(t):  Current  position  of  particle  i  at  iteration  t 

Pbestid-  Pbest  of  particle  i 

dbestid-  gbestof  particle i 

W:  Weight  function  for  velocity  of  particle  i 

C;:  Weight  coefficient  of  each  term 

The  first  term  on  the  RHS  of  the  equation  is  the  previous  velocity  of  the  particle.  Without  the  second  and  third 
term,  the  agent  will  keep  on  flying  in  the  same  direction  until  it  hits  in  the  boundary.  In  the  end  the  particles  will  try  to 
converge  to  the  Pbestor  gbest. 

Vmax:  So  from  the  convergence  Vmaxis  obtained  Vmax  that  is  the  maximum  change  in  one  particle  can  take  during 
the  iteration,  which  if  too  high  the  particle  must  fly  past  good  solutions,  similarly  if  it  is  too  small  particles  may  not  explore 
sufficiently  beyond  the  local  solution. 

C:  the  constant  d  and  C2  represent  acceleration  constants  or  learning  factors,  usually  considered  equal  and  ranges 
from  [0,4]  whenCl=C2=2 

W:  To  control  the  impact  of  previous  history  of  velocities  on  the  current  velocity  the  inertia  weight  W  is 
employed  which  influences  the  trade-off  between  global  and  local  exploration  abilities  of  the  flying  point.  Thus  suitable 
selection  of  inertia  weight  W  can  provide  a  balance  between  global  and  local  exploration  abilities  and  will  require  less 
iteration  on  average  to  find  the  optimum.  W  can  be  calculated  according  to  the  following  equation, 

W  =  Wmax  -  w^~w^  x  iter  (6) 

Itermax  isthemaximumnumberofiteratios 
Iter-is  the  number  of  iterations  until  the  current  stage 
Where  Wmax=  0.9  and  Wmin=  0.4 

Using  the  above  equation,  a  certain  velocity  that  gradually  gets  close  to  Pbest  and  gbest  can  be  calculated,  which  can 
be  further  used  to  calculate  the  current  position  with  the  following  equation 

Xid(t  +  l)=Xid(t)  +  Vid(t  +  l)  (7) 

Xid(t):  Current  position  of  particle  I  at  iteration  t 

Xid(t  +  1):  Current  position  of  particle  I  at  iteration  t+1 

Vid  (t  +  1):  Modified  velocity  of  particle  i 
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Algorithm 

The  following  PSO  algorithm  is  used  to  obtain  the  optimal  location  of  UPFC  for  loss  minimization. 

Step  1:  Initial  searching  points  and  velocities  are  randomly  generated  within  their  limits. 

Step  2:  Pbest  is  set  to  each  initial  searching  points.  The  best  evaluated  values  among  the  Pbest  are  set  to  gbest. 

Step  3:  New  velocities  are  calculated  using  equation  (5) 

Step  4:  lfVid(t+l)<Vdmin 

Then  Vid(t+l)  =  Vdminnnd  if  Vid(t  +  1)  > 

ThenVw(t  +  l)  =  Vdmax 

Step  5:  New  searching  points  are  calculated  using  equation  (7) 
Step  6:  Check  the  capacity  limits  constraints 
UPu(t  +  l)  >Pdmax 

Then  Pu(t  +  1)  =  Pdmax  and  if  Pu(t  +  1)  <  Pdmin 
TheaPid(t  +  l)=Pdmin 

Step  7:  Evaluate  the  fitness  value  for  new  searching  point.  If  evaluated  value  of  each  agent  is  better  than  previous 
Pbest  then  set  to  P^est -if  the  best  Pt,est  is  better  than  gbest\hen  set  to  gt,est 

Step  8:  If  the  maximum  iteration  is  reached  stop  the  process  otherwise  go  to  step  3. 


X  1 
pbest 


Figure  3:  Concept  of  Modification  Searching  Points  in  PSO 

RESULTS 

Simulation  studies  were  done  for  different  scenarios  in  IEEE  30  bus  power  system. 
Scenario  1:  Power  system  normal  operation  (without  FACTS  devices  installation). 
Scenario  2:  UPFC-GA  operation 
Scenario  3:  UPFC-PSO  operation 


The  first  scenario  is  normal  operation  of  network  without  installation  any  device.  In  second  UPFC-GA  was 
.  And  finally  UPFC-PSO  was  installed.  By  con 
improvement  and  reduces  the  losses  comparing  UPFC-GA. 


installed.  And  finally  UPFC-PSO  was  installed.  By  comparing  Scenario  2  &  3,  3ld  one  UPFC-PSO  was  given  best 
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Table  1:  Power  losses  with  and  without  UPFC 


SCENARIOS 

Total  Power 
Losses 

WithoutUPFC 

29.392 

UPFC  GA 

23.0900 

UPFC  PSO 

22.8591 
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Figure  4:  Voltage  Profile  with  and  without  Placement  UPFC 

CONCLUSIONS 

In  this  paper,  performance  of  transmission  system  with  UPFC  device  is  studied  using  Particle  Swarm 
Optimization  (PSO)  and  Genetic  Algorithm  (GA).  It  is  concluded  that  with  PSO  the  performance  studies  in  terms  of 
voltage  profile  and  power  loss  are  more  accurate  compared  with  GA.  Using  PSO  accurate  results  can  be  obtained.  Case 
Study  was  conducted  on  IEEE-30  bus  system. 
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